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PREFACE 

This  report  is  part  of  RAND's  continuing  efforts  on  the  detection 
of  high-altitude  nuclear  explosions .  It  Is  concerned  with  the  iono- 
spheric  propagation  of  hydromagnet ic  waves,  either  natural  or  resulting 
from  a  nuclear  burst.  The  results  of  previous  propagation  studies  are 
substantially  generalized  and  extended.  The  report  should  be  of  interest 
to  people  working  in  the  following  fields: 

(a)  Geomagnetism  and  related  Ionospheric  phenomena; 

(b)  Geomagnetic  effects  of  high-altitude  nuclear  bursts; 

(c)  Electromagnetic  wave  propagation  in  plasmas,  particularly  the 
ionosphere . 

The  author  is  a  consultant  to  The  RAND  Corporation. 
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In  this  paper  we  consider  the  two  cases  of  a  "fast”  or  "slow" 
hydromagnetic  wave  of  unit  amplitude  and  of  extremely  low  frequency 
(below  a  few  cps)  incident  on  the  upper  ionosphere  (below  the 
electron-density  maximum)  at  middle  or  high  latitudes,  and  we  study  its 
vertical  propagation  through  the  ionosphere.  From  approximate  solu¬ 
tions  of  the  propagation  equations ,  we  obtain  analytic  expressions  for 
l  magnetic  field  at  the  ground,  as  a  function  of  the  various  relevant 
quantities  such  as  wave  frequency,  polarization  of  the  incident  wave, 
ionospheric  ion  density,  ion  cyclotron  frequency,  ion-neutral  collision 
frequency,  and  the  dip  angle  of  the  earth's  magnetic  field.  The  results 
are  compared  with  those  of  previous  theoretical  studies  for  a  few  spe¬ 
cial  cases,  and  with  recently  reported  experimental  observations.  In 
both  cases,  the  agreement  is  good.  A  few  qualitative  features  of  our 
results  for  frequencies  between  about  .25  and  2  cps  are: 

(1)  For  middle  to  high  latitudes,  excluding  the  polar  regions, 
the  order  of  magnitude  of  the  transmitted  magnetic  field  is  insensitive 
to  the  polarization  of  the  incident  wave. 

(2)  The  daytime  transmitted  fields  should  be  predominantly  of 
right-handed  circular  polarization  (in  the  horizontal  plane)  regardless 
of  the  polarization  of  the  incident  fields. 

(3)  The  relative  daytime  amplitudes  of  the  transmitted  fields 
for  the  different  frequencies  are  insensitive  to  changes  in  the  iono¬ 
spheric  ion  density,  while  the  absolute  values  decrease  exponentially 
with  increasing  ion  density. 

(4)  Day-night  amplitude  ratios  vary  between  about  .05  and  1, 
depending  mostly  upon  daytime  ion  density,  and  are  insensitive  to  the 
wave  frequency  for  frequencies  above  about  .5  cps. 
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I .  INTRODUCTION 


In  this  paper  we  will  consider  theoretically  the  transmission  of 
extremely  low- frequency  (below  a  few  cps)  geomagnetic  pulsations  propa¬ 
gating  vertically  through  the  daytime  ionosphere  at  raddle  and  high  lat¬ 
itudes.  By  the  term  "ionosphere"  we  will  be  referring  to  the  region  be¬ 
low  the  F2  electron  density  maximum;  the  filtering  action  of  the  lower 
exosphere  will  not  be  considered  here. 

It  is  generally  assumed  that  most  of  the  extremely  low-frequency 
natural  fluctuations  in  the  geomagnetic  field  observed  at  the  earth's 
surface  have  their  source  somewhere  above  the  ionosphere  and  propagate 
downward  through  the  ionosphere  as  hydromagnetic  "waves."  If  we  are  to 
use  the  ground  data  to  make  deductions  concerning  the  .  mplitude,  fre¬ 
quency  of  occurrence,  and  frequency  spectra  of  such  disturbances  above 
the  ionosphere,  we  need  to  know  the  transmission  properties  of  the  iono¬ 
sphere  Itself.  Since  the  various  parameters  characterizing  the  iono¬ 
sphere  vary  with  latitude  and  fluctuate  with  season,  time  of  day,  and 
sunspot  conditions,  it  is  virtually  impossible  to  interpret  ground  data 
without  at  leas*  some  knowledgj  of  the  qualitative  dependence  of  the 
transmitted  fields  on  the  relevant  ionospheric  parameters,  and  of  the 
order  of  magnitude  of  the  changes  to  be  expected  in  the  transmitted 
fields  as  these  parameters  undergo  certain  changes. 

For  thir  purpose  it  is  very  useful  to  obtain  analytic  approxima¬ 
tions  of  the  transmitted  amplitudes  as  a  function  of  the  various  rele¬ 
vant  quantities  such  as  wave  frequency,  polarization  of  the  incident 
wave,  ionospheric  ion  density,  ion  cyclotron  frequency,  ion- neutral 
collision  frequency,  and  the  d*p  angle  of  the  earth's  magnetic  field. 

The  transmission  properties  of  the  ionosphere  for  hydromagnetic 
disturbances  in  the  frequency  range  of  interest  here  have  been  calcu¬ 
lated  for  a  few  special  cases.  Greifinger  and  Greifinger  (1965),  in  a 
paper  henceforth  referred  to  as  I,  obtained  an  analytic  expression  for 
the  ionospheric  transmission  coefficient  for  polar  propagation  through 
a  simple,  idealized  model  ionosphere.  Francis  and  Karplus  (1960),  and 
Karplus,  Francis,  and  Dragt  (1962),  in  two  works  henceforth  referred 
to  as  II  and  III,  respectively,  calculated  the  transmitted  amplitudes 
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of  hydromagnetic  waves  vertically  incident  on  the  ionosphere  at  an  al¬ 
titude  of  550  km  for  a  60-deg  dip  angle  of  the  earth's  magnetic  field. 
The  -esults  of  II  and  III  were  obtained  by  numerical  integration  of  the 
relevant  differential  equations  for  a  few  frequencies  in  the  range  of 
1  to  30  rad/sec  for  four  sets  of  ionospheric  parameters  corresponding 
tc  day  and  night  at  minimum  and  maximum  sunspot  conditions.  Wentworth 
(1964)  used  the  results  of  III  tc  obtain  an  empirical  expression  re¬ 
lating  attenuation  in  the  ionosphere  to  the  wave  frequency  and  to  the 
F9  peak-ion  density.  Prince  and  Bostick  (1S64)  and  Field  and  Greifinger 
(1965)  have  treated  the  case  of  vertical  propagation  at  the  equator. 

We  will  not  efer  again  to  these  last  two  works,  since  the  transmission 
characteristics  of  the  equatorial  ionosphere  are  quite  different  from 
those  at  other  latitudes,  and  we  are  concerned  here  only  with  middle- 
and  high- latitude  transmission. 

The  main  qualitative  features  of  the  results  of  I  for  daytime  polar 
propagation  are*. 

(1)  the  critical  dependence  of  ‘‘he  transmitted  amplitude  on  the 
polarization  of  the  incident  wave,  with  waves  of  left-handed  circular 
polarization  being  evanescent  in  the  lower  ionosphere  and  strongly  at¬ 
tenuated  for  frequencies  above  .1  to  #5  cps  (depending  upon  sunspot 
conditions),  and  waves  of  right-handed  circular  polarization  being  only 
■lightly  attenuated  for  frequencies  less  than  a  cps  (the  attenuation  in 
this  case  being  due  to  collisions!  absorption),  and 

(2)  the  existence  of  a  very  low-frequency  transmission  maximum  for 
a  wave  of  right-handed  circular  polarization.  (Unfortunately,  in  I,  the 
wave  modes  are  labeled  contrary  to  the  usual  convention.  The  electric 
(or  magnetic)  vector  for  the  wav  '  labeled  L  in  I  actually  has  right- 
handed  circular  polarization  when  viewed  along  the  earth's  magnetic 
field,  and  the  R  wave  of  I  has  left-handed  circular  polarization.) 

The  main  features  of  the  results  of  II  and  III  for  daytime  propa¬ 
gation  are: 

(1)  the  strong  dependence  of  the  transmitted  amplitudes  on  sunspot 
conditions, 

(2)  the  decoupling  of  the  fields  into  circularly  polarized  fields 
in  the  lower  ionosphere,  with  the  preferential  transmission  of  waves  of 
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right- handed  circular  polarization,  regardler.a  jf  the  polarization  of 
the  wave  incident  on  the  upper  ionosphere,  and 

(3)  the  insensitivity  of  the  order  of  magnitude  of  the  amplitude 
of  the  transmitted  magnetic  field  to  the  polarization  of  the  incident 
wave. 

For  the  case  of  nighttime  transmission  at  frequencies  below  a  few 
cps,  the  results  of  I,  II,  and  III  show  that  the  normal  ionosphere  is 
essentially  transparent  below  about  300  km;  however,  the  transmitted 
amplitudes  may  be  affected  by  anomalies  (such  as  sporadic  e-layer) . 
Reflection  from  the  earth,  of  course,  affects  nighttime  transmission 
as  it  does  the  daytime. 

in  this  paper  we  extend  the  enalyses  of  I,  II,  and  III  to  obtain 
approximate  analytic  expressions  for  the  transmitted  field  amplitudes 
of  low-frequency  hydromagnetic  waves  vertically  incident  on  the  daytime 
ionosphere  as  a  function  of  the  relevant  parameters.  We  give  numerical 
results  for  middle- latitude  parameters  obtained  from  the  ionosphere  pro¬ 
files  of  Sims  and  Bostick  (1963)  under  minimum  and  maximal  sunspot  con¬ 
ditions  and  predict  the  changes  if  different  values  are  employed.  More 
specifically,  we  obtain  approximate  rxpresriens  for  the  amplitude  of  the 
transmitted  daytime  magnetic  field  at  middle  and  high  latitudes  and  for 
its  ratio  to  the  amplitude  of  the  transmitted  nighttime  field  as  a  func¬ 
tion  of  wave  frequency,  polarization  of  the  incident  fields,  and  iono¬ 
spheric  ion  density.  In  addition,  we  investigate  the  effects  on  the 
transadtted  fields  of  the  altitude-dependence  of  the  ion— neutral  col¬ 
lision  frequency,  since  this  is  the  only  relevant  ionospheric  parameter 
which  changes  drastically  from  the  bottom  to  the  top  of  the  ionosphere. 

Our  results  are  compared  with  the  numerical  calculations  of  III 
and  with  experimental  observations  recently  reported  by  Campbell  (1967); 
in  both  cases  the  agreement  is  generally  good. 

The  main  qualitative  features  of  our  results  for  vertical  propaga¬ 
tion  of  frequencies  between  about  .25  and  2  cps  at  middle  latitudes  are: 

(1)  As  in  II  and  III,  the  order  of  magnitude  of  the  transmitted 
magnetic  field  is  insensitive  tc  the  polarization  of  the  incident  wave; 
changes  of  the  latter  are  not  likely  to  affect  the  transmitted  field  by 
more  than  a  factor  of  about  2. 
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(2)  The  amplitude  of  the  transmitted  field  decreases  exponentially 
vlth  increasing  density  of  ionospheric  ions — by  1/e  for  an  increased 
ion  density  of  about  55  percent. 

(3)  The  relative  daytime  transmitted  amplitudes  of  the  different 
frequencies  are  insensitive  to  changes  in  ion  density.  Thus  the  shape 
of  the  amplltude--frequency  curve  observed  at  the  ground  should  be  in- 
dsr indent  of  season,  time  of  day,  and  sunspot  conditions,  provided  that 
the  shape  of  the  amplitude — frequency  curve  in  the  upper  ionosphere  is 
unchanged .. 

(4)  The  ratio  of  the  day--night  amplitude  varies  between  .05  and  1, 
depending  upon  ion  density.  The  frequency-dependence  of  this  ratio  is 
slight  for  frequencies  above  about  .5  ops. 

(5)  The  fields  trensmitted  in  the  daytime  should  be  predominantly 
of  rlght-hcnded  circular  polarization  (in  the  horizontal  plane),  regard¬ 
less  of  the  polarization  of  the  incident  fields. 

If  for  polar  propagation  we  use,  instead  of  the  hort-range  expo¬ 
nential  form  employed  in  I,  a  more  realistic,  longer- rr.nge  function  for 
the  dependence  of  the  ion-neutral  collision  frequency  on  altitude,  the 
transmitted  amplitudes  become  smaller  than  those  of  I,  and  the  trans¬ 
mission  maximum  becomes  both  broader  and  flatter  than  that  of  I.  Other¬ 
wise,  our  results  for  polar  propagation  are  qualitatively  similar  to 
those  of  I. 
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II.  PROPAGATION  EQUATIONS 

For  the  aake  of  completeness,  we  shall  briefly  review  the  equations 
derived  in  II  governing  the  propagation  of  hydromagnetic  waves  in  the 
lo.tosphere.  Our  basic  assumptions  are  the  same  as  those  in  II;  the  iono- 
jhere  is  motionless  and  flat,  hydrostatic  and  gravitational  forces  may 
be  neglected,  the  motion  of  neutral  particles  and  the  collisions  between 
electrons  and  ions  may  be  neglected,  the  earth's  magnetic  field  is  uni¬ 
form,  the  parameters  describing  the  medium  vary  only  in  the  vertical 
direction,  and  the  normal  to  the  wave  is  in  the  verticrl  direction. 

We  also  linearize  the  equations  of  motion,  as  is  customary. 

We  shall  take  the  positive  z-directlon  to  be  vertically  downward 
and  the  static  uagnetic  field  to  be  in  the  x-z  plane.  The  electric 
field  E  and  thf  perturbed  magnetic  field  B  may  be  written 

E  -  (Ex,Ey,E8)  (1) 

B  -  <Bx,By,0) 

Assuming  a  time- dependence  of  the  form  e"iU)t,  and  neglecting  the  dis¬ 
placement  current,  we  st 'rt  with  the  wave  equation 

“*  s  -»  •+  2-*  Anita  ■— 

7  x  (7  X  E)  -  7(7  •  E)  -  7  E  -  -*-y=  J  (2) 

c 


where  J  is  the  current  density  and  is  related  to  E  by 


J  - 


o  e  (I 
o  o 


eo)  +  o 


itf- 


e  (E 
ov 


V3 


+  a,!  x  e 
2  o 


(3) 


The  vector  e  is  a  unit  vector  in  the  direction  of  B  ;  (e  )  ■  cos  <p, 

o  O  OZT’ 

and  (eD)x  “  sin  <p,  where  cp  is  the  complement  of  the  dip  angle.  The  con¬ 
ductivity  components  oq,  a^,  and  are  given  by 
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o  -  e  N. 

o  i 


t  m^Cv^  -  itu)  me('Je  "  i0J) 


.  A  | - +  — — j— r  I  <5> 

KC(vi  "  iW)2  +  V  °e[(ve  ■  iW)  +  “e]  ) 


v  -  itu 
e 


s.  t  “i 

N.\ —  2  2P 

(^C <V±  -  itu)  +  «!_] 


®  [(v  -  itu)2  +  tu2] 


whcie  N.  is  the  volume  density  of  ion  pairs,  and  mg  are  the  ion  and 
electron  masses,  and  are  the  ion  and  electron  collision  frequen¬ 
cies,  end  and  <ue  are  the  cyclotron  frequencies  of  ions  and  electrons, 
respectively.  As  in  II,  it  is  convenient  to  introduce  the  auxiliary 
quantities 


a,  -  a.  +  — 
3  1  (?, 


2  i  2 

O  -  cos  cp  +  —  sin  cp 

1  ao 


2  °3  2 

a/  -  cos  cp  +  —  sin  cp 

3  ao 


The  field  equations  for  the  transverse  components  Ex  and  Ey  of  the 
electric  field  are  then  given  by 

d2  .  E  +  4nl-  °2  C°S  --  E  -  0  (10) 

t?+  c2  °1#JX  c2  y 


d2  4TTitu  al  /  E  4nitu  °2  C°J_*  .  0 
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The  longitudinal  component  E  of  the  electric  field,  which  is  not  a 

£ 

dynamic  variable,  can  be  expressed  in  terms  of  and  E^: 

E  -  (e  —  -f  E  cos  q>(-i  -  l)l 

z  ai  1  y  °o  x  Voo  n 


Henceforth,  we  shall  concern  ourselves  only  with  field  components  in 
the  x-y  plane. 

An  examination  of  representative  height  profiles  of  collision  and 
cyclotron  frequencies  shows  that,  except  for  very  small  dip  angles  (a 
few  degrees  or  less),  electron  collisions  are  negligible  and 


c '  ~  ~  cos  '  «p 


Substituting  Eq.  (13)  into  Eqs.  (10)  and  (11),  and  neglecting  ve> 


ve  get 


d2  .  4rriu>  ffl 

— 2 +  2 — r~ 

.dr  c  coa '  cp 


e  +»^s-!2_e 

X  2  cos  cp  y 


d  ,  4ni<»  _  „  4TTiu>  °2 

— x  +  — =—  a.  E - - - E 

da  c2  1J  y  c2  co8  *  x 


where 


4TTit«  (uM  u)(u?  +  ivj) 


2  2  2 
va  +  (vt  -  iw) 


4TTitU  0, 


.  /M  (vi  • lal) 

v?  \V  I  (vt  -  i»)2  + 
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V 

a 


4rNiMi 


1/2 


-  Alfven  speed 


(18) 


One  boundary  condition  is  that  the  field  components  E  and  E  van- 

x  y 

lsh  at  the  surface  of  the  earth,  which  we  treat  as  a  perfect  conductor 
at  the  low  frequencies  of  interest  here.  Additional  boundary  conditions 
will  be  discussed  later. 
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III.  IONOSPHERIC  REGIONS  AND  NORMAL  MODES  OF  RADIATION 


For  hydromagnetic  waves  in  the  frequency  regime  below  few  cps, 
the  ionosphere  consists  of  three  regions  that  are  characterized  by  the 
relative  magnitudes  of  v^»  uj^,  and  u),  and  by  cos  cp  (see  I).  In  order 
of  increasing  altitude,  these  regions  are: 

(a)  Hall  region  (a>  <  cos  cp) 

(b)  Transition  region  (uo  <  <  tu^/cos  cp) 

(c)  Alfven  region  (v^  <  <  (u^) 

The  location  if  the  "boundary"  between  the  Alfven  and  transition 
regions  depends  upon  <u  and  upon  sunspot  conditions.  For  a  given  <u,  the 
altitude  of  this  boundary  increases  with  the  sunspot  number  as  the  re¬ 
sult  of  the  higher  collision  frequencies  associated  with  periods  of  high 
solar  activity. 

The  height  of  the  upper  "boundary"  of  the  Hall  region  is  Indepen¬ 
dent  of  <u  and  only  slightly  dependent  upon  latitude  at  high  and  middle 
latitudes.  (At  very  low  magnetic  latitudes,  where  cos  9  is  saall,  the 
Hall  region  becomes  thin,  and  it  vanishes  near  the  equator.)  The  thick¬ 
ness  of  the  Hi'.ll  region  also  depends  upon  sunspot  conditions,  the  upper 
boundary  increasing  in  altitude  from  about  120  km  at  sunspot  minimum  to 
about  130  km  at  sunspot  maximum. 

The  character  of  the  propagation  equations  and  the  properties  of 
the  normal  modes  of  the  radiation  are  quite  different  in  the  different 
regions.  In  the  Hall  region,  |aJ/cos  9<  |a2J.  If  loj/cos  9  is  negli¬ 
gible  compared  to  [cjjl,  Eqs.  (14)  to  (18)  tecome  approximately 


dz 


21Ex  *  iEy 


F  —  [E  T  iE  ] 
uij  v‘  x  yJ 

a 


(19) 


The  two  normal  modes  (E  ■  ±  iE  )  are 

y  x 

horizontal  plane,  and  are  characterized  by 


circularly  polarized  in  the 
distinct  indices  of  refraction. 
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i 
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one  real  and  the  other  imaginary,  within  the  approximation  used.  Thus, 
the  two  modes  are  sometimes  referred  to  as  propagating  (p)  and  evanes¬ 
cent  (e).  When  viewed  along  the  earth's  magnetic  field,  the  E  (or  B) 
vector  rotates  clockwise  in  the  p-mode  and  counterclockwise  in  the  e- 
mode.  In  I  the  p-  and  e-modes  are  unfortunately  referred  to  as  L  and 
R  modes,  respectively,  contrary  to  the  usual  convention.  This  termi¬ 
nological  confusion  results  from  the  dependence  of  the  sense  of  rota¬ 
tion  on  the  direction  from  which  the  fields  are  viewed,  and  we  shall 
henceforth  use  the  unambiguous  "p"  and  "e"  labels  for  the  modes. 

In  the  Alfven  region,  Eqs.  (14)  to  (18)  become 


When  V  and  u>.  are  constant,  solutionr  have  the  form  e  ,  and 
a  1 

the  dispersion  relationsnip  is  given  by 


where 


« •  8 


1  + 


2 

cos  9 


2 

cos  9 


4 


2  2 

.4  .  4u>  cos  <p 

sin  9  + - - 1 


(22) 


(23) 


is  the  phase  velocity. 
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III.  IONOSPHERIC  REGIONS  AND  NORMAL  MODES  OF  RADIATION 


For  hydromagnetic  waves  in  the  frequency  regime  below  a  few  cps, 
the  ionosphere  consists  of  three  regions  that  arc  characterized  by  the 
relative  magnitudes  of  and  u>,  and  by  cos  cp  (see  I).  In  order 

of  increasing  altitude,  these  regions  are: 

(a)  Hall  region  (cu  <  cos  cp) 

(b)  Transition  region  (cw  <  <  tu^/cos  cp) 

(c)  Alfven  region  ^  <  cu  <  cu^) 


The  location  of  the  "boundary"  between  the  Alfven  and  transition 
regions  depends  upon  cu  and  upon  sunspot  conditions.  For  a  given  cu,  the 
altitude  of  this  boundary  increases  with  the  sunspot  number  as  the  re¬ 
sult  of  the  higher  collision  frequencies  associated  with  periods  of  high 
solar  activity. 

The  height  of  the  upper  "boundary"  of  the  Hall  region  is  indepen¬ 
dent  of  cu  and  only  slightly  dependent  upon  latitude  at  high  and  middle 
latitudes.  (At  very  low  magnetic  latitudes,  where  cos  cp  is  small,  the 
Hall  region  becomes  thin,  and  it  vanishes  near  the  equator.)  The  thick¬ 
ness  of  the  Hall  region  also  depends  upon  sunspot  conditions,  the  upper 
boundary  increasing  in  altitude  from  about  12C  km  at  sunspot  minimum  to 
about  130  km  at  sunspot  maximum. 

The  character  of  the  propagation  equations  and  the  properties  of 
the  normal  modes  of  the  radiation  are  quite  different  in  the  different 
regions.  In  the  Hall  region,  jajj/cos  cp  <  |a2|  .  If  [aj/cos  cp  is  negli¬ 
gible  compared  to  | CT2 1  »  £<18-  (14)  to  (18)  become  approximately 


A 


E  ?iE 
x  y 


*  —  -4  [E  T-  iE  ] 

»i  v2  *  y 

a 


(19) 


The  two  normal  modes 
horizontal  plane,  and  are 


(E  =  ±  iE  )  are  circularly  polarized  in  the 
y  x 

characterized  by  distinct  indices  of  refraction. 
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one  real  and  the  other  imaginary,  within  the  approximation  used.  Thus, 
the  two  modes  are  sometimes  referred  to  as  propagating  (p)  and  evanes¬ 
cent  (e) .  When  viewed  along  the  earth's  magnetic  field,  the  E  (or  B) 
vector  rotates  clockwise  in  the  p-mode  and  counterclockwise  in  the  e- 
mode.  In  1  the  p-  and  e-modes  are  unfortunately  referred  to  as  L  and 
R  modes,  respectively,  contrary  to  the  usual  convention.  This  termi¬ 
nological  confusion  results  from  the  dependence  of  the  sense  of  rota¬ 
tion  on  the  direction  from  which  the  fields  are  viewed,  a wi  we  *hall 
henceforth  use  the  unambiguous  **pM  and  "e"  labels  for  the  modes. 

In  the  Alfven  region,  Eqs.  (7.4)  to  (18)  become 


When  V  and  u).  are  constant,  solutions  have  the  font  e  ,  and 
*■  i 

the  dispersion  relatlr  nship  is  given  by 


V2k2  V2 

~V  s~r 

■  v 


where 


cos  Cp 


sin  cp  + 


4a?  cos  cp  I 


(22) 


(23) 


is  the  phase  velocity. 
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Thus,  Eqs.  (30)  and  (31)  provide  a  smooth  transition  from  linearly 
polarized  to  circularly  polarized  normal  modes  as  cp  decreases  from 
about  y/ 2uo/ cu^  to  0. 

The  nomenclature  of  the  two  modes  in  the  Alfven  region  is  •_  ource 
of  as  much  confusion  as  th  nomenclature  of  the  modes  in  the  Hall  re¬ 
gion.  In  II  and  III,  the  fast  and  slow  waves  are  referred  to  as  "ordi¬ 
nary"  and  "extraordinary,"  respectively.  Other  authors  reverse  this 
nomenclature,  referring  to  the  fast  wave  as  "extraordinary."  Other 
terms  commonly  used  for  the  slow  wave  are  "Alfven  mode,"  "L-mode,"  and 
"ion  cyclotron  wave,"  (the  last  two  referring  to  the  wave  properties 
for  parallel  propagation).  The  fast  wave  is  sometimes  described  as 
"R-mode"  or  "isotropic;"  the  latter  term  refers  to  the  lack  of  depen¬ 
dence  of  its  phase  velocity  on  cp,  and  the  former  to  its  polarization 
for  parallel  propagation.  Heedless  to  say,  it  would  be  helpful  if  a 
standard  nomenclature  were  adopted. 

We  have  seen  that  the  waves  can  be  described  in  terms  of  (nearly) 
uncoupled  normal  inodes  in  idealized  Alfven  and  Hall  regions.  The  po¬ 
larization  of  the  two  modes  in  the  Alfven  region  is  elliptical,  the  two 
ellipses  becoming  mutually  perpendicular  straight  lines  when  cp  »  V Zm/cu^ 
and  circles  when  cp  «  v^uu/u^.  In  the  Hall  region,  the  normal  inodes  are 
circularly  polarized  in  the  horizontal  plane. 

The  transition  from  one  pair  of  normal  modes  in  the  Alfven  region 
to  another  pair  in  the  Hall  region  takes  place  in  the  transition  region, 
where  the  waves  are  coupled  because  of  collisions.  This  coupling  of 
modes  is  an  important  property  of  the  daytime  ionosphere,  and  leads  to 
changes  in  the  polarization  of  the  radiation  as  it  passes  through  the 
ionosphere. 

An  exception  occurs  when  the  wave  propagates  parallel  to  the  mag¬ 
netic  field.  For  cp  >»  0,  and  for  a  model  ionosphere  in  which  all  parame¬ 
ters  vary  only  in  the  vertical  direction,  the  two  circularly  polarized 
normal  modes  remain  uncoupled  as  the  radiation  descends  vertically  from 
the  Alfven  region  to  the  Hall  region.  The  fast  mode  of  the  Alfven  re¬ 
gion  becomes  the  propagating  mode  of  the  Hall  region,  and  the  slow  mode 
of  the  Alfven  region  becomes  the  evanescent  mode  of  the  Hall  region. 

This  was  the  case  in  I,  where  polar  propagation  (cp  =  0)  was  considered. 
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When  <p  »  V 2uj/oj^ ,  the  radiation  in  the  Hall  region  consiats  of  a 
mixture  of  p  and  e  modes,  regardless  of  the  polarization  of  the  inci¬ 
dent  wave.  The  extent  to  which  the  transmitted  fields  depend  upon  the 
wave  mode  of  the  incident  radiation  will  be  discussed  later. 
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IV.  PHYSICAL  MODEL 


In  order  to  find  analytic  solutions  of  the  wave  equation,  we  need 

to  establish  the  functional  dependence  of  the  parameters  tu, ,  V  ,  and 

X  b 

on  z. 

According  to  the  profiles  of  Sims  and  Bostick  (1963),  varies 
little  with  altitude  in  the  ionosphere,  although  it  is  latitude-de¬ 
pendent.  For  each  latitude  considered,  we  will  assume  <jj^  cO  be  inde¬ 
pendent  of  z. 

The  Alfven  speed  Vfl  depends  upon  latitude,  time  of  day,  sunspot 
conditions,  and  altitude.  When  considered  as  a  function  of  'ncreaslng 
altitude  from  the  surface  of  the  earth,  Vfl  decreases  abruptly  at  an 
altitude  of  about  100  km.  We  refer  to  100  km  as  the  "bottom"  of  the 
ionosphere,  since  the  region  below  100  km  is  characterised  by  values 
of  Vfl  so  large  that  the  medium  behaves  like  a  vacuum  for  the  frequen¬ 
cies  of  interest  here.  For  altitudes  above  100  km  but  below  the 
minimum  of  the  Alfven  speed,  V  decreases  with  altitude,  the  total  de- 

fl 

crease  from  the  bottom  to  the  top  of  the  ioncsphere  being  about  a  fac¬ 
tor  of  2  or  3,  depending  upon  sunspot  conditions.  In  the  Hall  region 
the  variations  are  slight,  and  V  can  be  taken  as  constant.  Similarly, 

B 

in  the  Alfven  region  decreases  slowly  up  to  the  F2  minimum  and  can 

be  regarded  as  constant.  In  the  transition  region,  we  shall  see  that 

2 

V  enters  the  equations  only  as  part  of  the  factor  v./V  ;  thus  we  will 

fl  X  B 

not  need  to  be  concerned  with  V  itself. 

a 

The  ion-neutral  collision  frequency  changes  by  a  few  orders 

of  magnitude  from  the  lower  to  the  upper  portions  of  the  ionosphere, 

decreasing  from  several  thousand  collisions/sec  at  100  km  to  a  value 

between  about  .1  and  1  collision/sec,  depending  upon  sunspot  conditions, 

at  about  300  km.  We  shall  be  concerned  with  the  functional  dependence 

2 

on  z  of  the  quantity  v4/V  in  the  transition  region;  this  quantity  is 
usually  expressed  in  terms  of  a  scale  height  Q: 


a 


(32) 
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'•here  we  have  taken  z  -  0  to  be  the  altitude  for  which  -  ou^.  The 
!e  height  £  Increases  with  altitude  from  the  bottom  to  the  top  of 
lw-.:'ephere,  varying  from  a  few  km  near  the  bottom  to  about  100  km  near 
the  top.  Since  the  height  of  the  upper  boundary  of  the  transition  re¬ 
gion  Is  frequency-dependent,  extending  to  nearly  t^e  top  of  the  Iono¬ 
sphere  for  the  lowest  frequencies  considered,  we  need  an  approximation 
of  £(z)  that  Is  applicable  to  the  entire  Ionosphere  above  the  Hall  re¬ 
gion. 

In  I,  £  was  assumed  to  be  constant  and  equal  to  Its  value  £q  at 
z  ■  0.  This,  of  course,  Is  the  simplest  approximation,  but  since  It 
greatly  underestimates  the  collision  frequencies  at  high  altitudes,  It 
reads  to  an  overestimate  of  the  amplitude  of  the  transmitted  field  for 
those  frequencies  at  which  the  transition  region  extends  a  distance 
greater  than  about  2 £q  above  z  ■  0.  For  frequencies  high  enough  so 
that  the  thickness  of  the  transition  region  Is  of  order  2£q  (or  less), 
£  can  be  replaced  by  £q. 

The  simplest  way  to  approximate  £  to  take  account  of  Its  Increase 

with  altitude  Is  to  assume  It  to  be  a  linear  function  of  z.  This  leads 

to  a  power  law,  rather  than  an  exponentl  '.1  dependence  of  v  /V2  on  al- 

1  8 

tltude: 


a 


UL' 


<va>  z-0 


(33) 


Equation  (33)  provides  a  reasonably  good  approximation  to  v  Vfl  when 
the  power  P  Is  equal  to  about  2  and  £q  Is  about  10  km  at  sunspot  mini¬ 
mum  and  15  km  at  sunspot  maximum.  This  can  be  seen  from  Fig.  1,  where 

2 

we  have  plotted  Eq.  (33)  together  with  the  values  of  v. /V  given  by 

i.  8 

Sims  and  Bostick  (1963).  The  latitude-dependence  of  the  functional 
2 

form  of  v./V  is  slight;  Eq.  (33)  is  equally  applicable  at  all  lati- 
tudes  considered  here.  We  shall  use  Eq.  (33)  with  P  =  2  and  compare 
the  results  with  those  obtained  by  using  a  constant  value  of  £,  as  In 
I. 
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-  fo  =  10km  at  sunspot  minimum 

f0  =  15km  at  sunspot  maximum 
-  -z=altitude -100-2  f0(km) 


i\  D-z/2rJ2 


Sunspot 

minimum. 


Values  obtained  from 
the  profiles  of  Sims 
and  Bostick  ( 1963) 


r  Sunspot 

jQ  maximum. 


5-/7; 


100  140  ISO  220  260  300  340 

Altitude  (km) 


Fig.  1  --  The  ion-neutral  collision  frequency  divided  by  the  square 
of  the  Alfvt*  speed,  normalise  1  to  unity  at  the  altitude  at  which 
the  ion-neutral  collision  f.^naency  equals  the  ion  cyclotron 
frequency  ui^.  The  data  points  are  from  the  ionospheric 
profiles  of  Sims  and  Bostick  (1963)  for  a  magnetic 
dip  angle  of  74°.  The  solid  curves  are  the 
empirical  fits  to  the  data  points  given 
by  Eq.  (33)  when  P  ■  2. 


aw#*#® 
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V.  BOUNDARY  CONDITIONS  AND  TRANSMISSION  MATRIX 


Before  we  find  solutions  to  the  propagation  equations,  we  must 
specify  the  boundary  conditions.  The  vanishing  of  the  electric  field 
at  the  surface  of  the  earth  gives  as  one  boundary  condition 


E  (z.) 

x,y  1 


(34) 


where  is  the  coordinate  of  the  bottom  of  the  ionosphere  and  h  is  the 
thit  kness  of  the  vacuum  region  (100  km) . 

The  choice  of  the  additional  boundary  conditions  is  rather  arbi¬ 
trary,  since  we  are  disregarding  the  source  of  the  radiation.  Most 
source  models  ;urrently  assume  that  the  radiation  impinges  on  the  iono¬ 
sphere  either  as  a  slow  wave  arriving  from  a  source  in  the  magnetosphere 
or  as  a  fast  wave  that  has  arrived  from  another  latitude  via  an  iono¬ 
spheric  duct  located  near  the  F2  electron-density  maximum.  Thus  for 
convenience  we  shall  consider  the  two  cases  where  the  incident  wave  is 
either  a  pure  slow  (S)  or  a  pure  fast  (F)  wave. 

In  order  to  translate  this  assumption  into  a  specific  boundary  con¬ 
dition  for  the  incident  fields,  let  us  examine  the  asymptotic  form  of 
the  electric  and  magnetic  field  components  in  the  Alfven  region.  In 


this  region  the  components  of  the  electric  field  satisfy  Eqs .  (20)  and 
(21).  We  will  assume  that  |u)/uj^|  «  1,  that  |(tu/u)i)ExJ  is  negligible 
compared  to  |e  |,  and  that  |  (uo/uj  ) E  |  is  negligible  compared  to  |e  |. 

y  i  y  x 

Equations  (20)  and  (21)  then  become  uncoupled  into  the  two  equations 


(35) 
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From  che  induction  equation,  we  know  that  B  satisfies  Eq.  (35)  and 
satitfies  Eq.  (35).  Therefore,  the  field  components  in  the  Alfven 
region  can  be  written  in  the  form 


ik  z  -ik  2 

B  -  AeX+(rA+rA)e  X 
x  x  xx  x  xy  yy 


where 


ik  2  -ik  2 

B  -Ae^+(rA+rA)e  ^ 

y  y  '  yy  y  yx  x' 


.  ik  2  -ik  2 

E  “  —  Ae  X  -  (r  A  +  r  A)e  x 
y  c[x  xx  x  xy  y 


E  ■  c.°3  SL  e  y  +(r  A  +  r  A)e 

x  c|y  yyy  vx  x 


-ik  2] 

y 


k  E  xT 
x  Va 


V  C08  CP 
a 


We  have  expressed  the  reflected  wave  in  terms  of  the  elements  r  , 

xx’ 

r  , . . . ,  of  a  2-by-2  matrix.  The  presence  of  nonvanishing  off-diagonal 
xy  , 

elements  r  and  r^  is  a  consequence  of  coupling  below  the  Alfven  re¬ 
gion.  A  particular  field  component  may  be  present  in  the  reflected  wave 
even  if  it  is  absent  from  the  incident  wave.  (Thus  Bx  contains  a  term 
proportional  to  A  ,  etc.) 

The  reflection  matrix  is  unknown  and  can  be  determined  only  ^ter 
Eqs .  (14)  and  (15)  are  solved.  In  order  to  apply  the  boundary  condi¬ 
tion  that  specifies  the  mode  of  the  incident  wave  in  the  Alfven  region, 
it  is  convenient  to  eliminate  the  unknown  reflection  terms  from  Eqs. 

(37)  to  (40)  and  to  relate  the  field  components  and  their  derivatives 
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to  specified  components  of  the  incident  field.  From  Eqs.  (37)  to  (40) 
we  get 


where 


dE 

E 


iS" 

V 


dz 


-  cos  cp(B  ). 
c  T  y'inc 


i  dE 
E 

y  ia)  dz 
V 

a 


"c  <Vlr.c 


(B  ),  ■  A  e 

'  x7 inc  x 


Ik  z 
x 


<Vinc 


=  A  e 

y 


ik  z 

y 


(43) 


(44) 


(45) 


(46) 


The  7  or  S  character  of  the  incident  wave  determines  the  polariza¬ 
tion  of  the  incident  fields  and  gives  the  relative  values  of  A  and  A 

°  x  y 

as  a  function  of  cp.  The  magnitudes  of  A  and  A  will  be  chosen  so  as 

x  y 

to  normalize  the  incident  magnetic  vector  to  unit  magnitude,  i.e., 

I*J2  +  I  Ay  1 2  -  i  <«> 

Thus,  combining  Eqs.  (24)  to  (31)  with  Eqs.  (39),  (40),  and  (47),  we 
obtain: 

Incident  Fast  Wave 


A 


y 


A 

x 


0, 


1 


for 


(48) 
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Incident  Slow  Wave 


Although  Eqs .  (50)  and  (53)  are  the  limits  of  Eqs.  (49)  and  (52),  re¬ 
spectively,  as  <p  -*  0,  it  is  useful  to  consider  parallel  propagation 
separately,  since  it  is  an  important  special  case  in  magnetohydrody- 
namic-wave  theory. 

Thus  Eqs.  (43)  to  (46),  in  conjunction  with  Eqs.  (48)  to  (53), 
constitute  our  boundary  conditions  in  the  Alfven  region. 

The  transmitted  fields  are  most  conveniently  described  in  terms 
of  circularly  polarized  fields,  regardless  of  the  polarization  of  the 
incident  fields,  because  of  the  effect  of  the  Hall  region  on  the  radi¬ 
ation.  Following  III,  we  define  a  magnetic  transmission  matrix  T  which 
relates  the  magnetic  field  at  the  ground  to  the  incident  field: 
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j  =  [B  -  IB  ]_  . 

Ground  x  yJGround 


VT  [Tx<w,cp)(Bx)inc  +  Ty(u>,cp)(By)lnc]  (54) 


For  positive  frequencies,  the  matrix  T  describes  the  transmitted 
p-mode,  and  for  negative  frequencies,  the  transmitted  e-mode.  The  in¬ 
dividual  field  components  E?  and  By  and  the  total  field  magnitude  |b| 
at  the  ground  are  related  to  ^-roun<j  By 


(2BJ 


x  Ground  Ground 


<2By) 


Ground 


'Ground 


Ground ' 


'Ground  ^  ^ 

(55) 

BGround(-“)] 

(56) 

IiHJLL 

■"Ground 

(57) 

Therefore,  if  T  (u»,q>)  and  T  (u>,9)  are  determined  '  ’  can  determine 
x  y 

the  transmitted  magnetic  field  for  any  A  and  A  ,  and  it  would  seem 

x  y 

unnecessary  to  apply  any  restrictions  such  as  Eqs.  (48)  to  (53)  to  the 

values  of  A  and  A  . 

*  y  _ 


However,  when  9  »  VZuu/u^,  we  will  see  that  the  methods  of  ap¬ 
proximation  to  be  used  in  solving  the  equations  will  rely  on  the  as¬ 
sumption  that  either  A  or  A  vanishes,  so  that  we  will  be  calculating 

x  y 

the  values  of  T  and  T  under  these  conditions.  The  transmitted  mug- 
x  y 

netic  field  under  other  conditions  of  polarization  of  the  incident 
radiation  can  then  only  be  deduced  by  means  of  an  "educated  guess." 

For  example,  if  the  Incident  radiation  is  unpolarized,  it  would  seem 
reasonable  to  assume  that  the  transmitted  magnetic  field  would  be  ap¬ 
proximately  equal  to  the  average  of  the  two  values  that  result  from 

the  condition  that  either  A  or  A  is  zero. 

x  y 

We  will  see  that  It  |  and  |t  |  will  fortunately  be  of  the  same 
x  y 

order  of  magnitude  for  middle  latitudes,  so  that  the  transmitted  field 
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magnitudes  will  not  be  particularly  sensitive  to  the  polarization  of 
the  incident  radiction. 

When  cp  is  small  (high  magnetic  latitudes),  the  propagation  equations 

become  independent  of  cp,  and  |t  |  and  |t  j  become  equal.  In  this  case, 

x  y 

we  can  indeed  find  the  transmitted  magnetic  field  for  arbitrary  A and 
Ay,  and  can  dispense  with  the  restriction  that  these  quantities  satisfy 
the  conditions  appropriate  to  an  Incident  F-  or  S-wave.  However,  since 
in  some  models  it  well  may  be  that  the  source  generates  waves  of  a 
given  mode,  such  a  boundary  condition  still  appears  sensible. 


-24- 


VI .  APPROX  ATIONS  AND  SOLUTIONS 


To  solve  our  equ&tions,  let  us  first  rewrite  Eqs.  (14)  to  (17)  in 
terms  of  dimensionless  variables  (denoted  by  a  tilde)  with  V  ,  and 
Vfl/u)i  at  z  »  0  as  the  units  of  frequency,  speed,  and  length,  rerpectively . 


u>  +  iuov . 


s~2 

dz 


E  - 

x 


— 2 
iuh^ 


9*  COB *  cp(l  +  ~j)J  x  9^  cos  <p(l  +  vt) 

8  X  8  X 


E  - 

y 


(58) 


.2  4-  i<uv. 

4.  _ _ i 

~2  «2  .  -2 


dz 


iuh;. 


E  + 

y 


9‘(1  +  vJ)J  y  9*  COS  vp(l  +  v*) 


(59) 


When  <p  <  V 25,  we  can  replace  cos  cp  by  1,  since  we  are  neglecting 
terms  of  order  5  compared  to  " ,  end  thus,  as  in  I,  Eqs.  (58)  and  (59) 
can  be  combined  into  one  equati  for  the  variable  (Ex  -  iE^) : 


where 


and 


0 


E  -  iE 
x  y 


(60) 


(61) 


e 


(62) 
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For  cp  >  V25*,  we  can  also  write  an  equation  for  £  in  either  of  two 


identical  forms 


4 + 4 

Idz  ^>y 


e  “  Q*,y(*> 


where 


~2  ~~ 
+  ilU'J. 


2/  1<3i 

( 1  -4-  - r  cos  cp 

*\  ?2 

'  a 


~2 

uj  +  iu>v. 


«2  2 
V  cos  cp 
a  Y 


\  V  cos  cp 


Qx(z) 


in  . 

T!  <V1  C0‘  ®  E„  -  Ex> 


(1  +  Vj/cos  cp 


i(l  +  cos 


Qy(«) 


i_8inii£/E  +  2!iJL_V 

1  +  V?  ‘  y  C0S  V 


^  cos^  cp[l  + 


„2 

V  cos 
a 


For  convenience  of  solution,  we  choose  the  subscript  x  or  y  in  Eq.  (63) 
according  to  the  boundary  conditions  on  the  incident  wave  (the  sub¬ 
script  x  for  A  -  1,  A  »  0,  and  y  for  A  ■  1,  A  =0).  _r>  the  limit 

x  y  y  x 

as  cp  -*  0,  «  -  q  and  Q  -•  0,  so  that  Eq.  (63)  reduces  to  Eq.  (60)  . 

Xj  y  o  x ,  y 

The  magnitude  of  the  inhomogeneous  term  Q  is  small  relative  to 

1  2  1  y 

|q  F|  in  the  Hill  region  and  in  the  upper  transition  and  lower  Alfven 

2 

regions.  The  magnitude  of  Q  relative  to  |q  t\  is  largest  near 
z  -  0(~  <*s  1),  where 
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x»y 


sin  epl  q_  £| 
Y  x,y  ' 


(68) 


At  middle  and  (especially)  at  high  latitudes,  sin  cp  is  small. 

(For  example,  at  an  intermediate  magnetic  latitude  of  about  40°, 

2 

sin  <p  «  .25).  Therefore  0  (z)  represents  a  small  perturbation  ex- 

j.  ,  y 

tending  over  a  small  region  near  z  «  0,  and  in  view  of  the  approximate 

2 

nature  of  the  analytic  representation  of  v. /V  which  we  shall  employ, 

1  3 

it  seems  pointless  to  retain  Q  (z)  in  Eq.  (63). 

—  x ,  y 

If  we  neglect  Q  (z),  Eq.  (63)  then  becomes  a  homogeneous  equa- 

x.y 

tion  of  the  same  form  as  Eq.  (60): 


d2£ 


dz 


~2  +  *£,3 


(69) 


To  solve  Eq.  (69)  it  is  necessary  to  have  an  analytic  approxima- 
~  «2  2 

tion  for  the  quantity  v,/v  contained  in  q”  .  To  study  the  sensitivity 

la  *>y 

of  the  results  to  different  assumptions,  we  shall  consider  both  a  short- 
range  exponential  dependence  on  altitude  (as  in  I) 


a 


(70) 


and  a  long-range,  inverse-square  dependence 


li 

v2 

a 


z  >  2 C  MZ. 
o  1 


(71) 


It  is  necessary  to  find  solutions  of  Eq.  (69)  for  the  range  of  alti¬ 
tudes  z  i  zo(tu)  Cwhere  *  ($)  lie**  in  the  lower  Alfven  region  near 

the  altitude  at  which  ~  U)]  over  which  collisions  are  important. 
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SOLUTIONS  FOR  VERY  SMALL  u 

In  the  limit  as  u)  -•  0,  £  becomes  a  linear  function  ox  z: 


£  - 


2k  h 
_ XjJL. 


1  + 


zi  - 2 


1  •  *  2!  -  2o> 


[Bx  -  lBy  cos  <p|lnc  (72) 


For  small  but  finite  3,  £  will  be  a  nearly  linear  function  of  z.  To  ob¬ 
tain  the  next  approximation,  let  us  write  Eq.  (69)  as  an  integral  equation: 


6(z)  - 


”  z  -  z 

2k  h 

i  +  -U— 

x.y 

h 

1  -  ll^x>y(h  +  \  -  *0) 

B  -  iB  cos  cp 
x  y  T 


inc 


+  G(z,z /)q^^,(z  ^^(z ')dz # 


(73) 


where  the  Green  function  G^.z*)  is  given  by 

hU^zjU^z') 


G(z,z  *)  - 


/  2n  -  *  \ 
1  -  ik  h  (1  +  — -  °  ) 

x»y  \  h  / 


hU2(z)U1(z/) 


/  zT  -  S  \ 

-  ik  h  (l  +  — -  °) 

x*y  v  h  / 


z  £  z  £  z ' 
o 


(74) 


where 


U.  (3)  -  1  + 


Z1  -  Z 


U2(z)  -  1  -  ik(z  -  lQ) 


(75) 


(76) 


uJi lii  k.i:\.L-.s.uniUu  la^^ii-iuiiiutukLitLli^i.uuvwii-ajUl-itlfihai'iittfiili 
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When  |u>|  is  so  small  that  j<uzo(ui)|  <  1,  the  integral  can  be  estimated 
easily  by  replac  ng  U2(i’;)  by  1  and  £(z/)  by  a  linear  function  of 


which  we  write  in  the  form 


M  / 

/  Z,  -  Z  , 


Solving  Eq.  (73)  for  £(z^) ,  we  obtain 


2k  h 

— [B  -  iB  cos  cp] 
c  x  y  TJinc 

IU  11  * 

-f1/  *],  "  *\  2  ~~ 

1  -  h  i  l1  + 


Th  integration  is  trivial  for  either  the  power- law  or  exponential  form 
of  v1/Vfl.  From  fi(z^)  we  can  obtain  (B)groun<j  by  applying  the  induction 
equation  at  the  earth: 


Ground 


c/de\ 

'  *«W„ 


Ground 


From  (B)Ground  an''  Eq.  (54,,  we  can  then  find  the  magnetic  transmission 

matrix  T,  the  results  for  which  are  as  follows  for  the  two  different 
2 

forms  of  v. /V  : 

i  a 


—  )(l  +  iy,) 

i  cp/  Tly 
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where 


The  magnetic 
limit  |  v^ujz^  | 
When 

(b) 


where 

Y2 


to 


C3 


TJ  fo 

'l  °  2z, 


(81) 


to  .7  for  sunspot  minimum 
ri  .8  for  sunspot  maximum 

transmissi  n  coefficient  of  I  reduces  to  Eq.  (80)  in  the 

«  1. 


(‘  - 1) 


T  - 

x 


_ VT _ 

2.5 

i - +  i- 

(cos  cp)t/Z 


(82) 


T  - 


_ -  y/2i _ 

2.5  t  Sfi 

1  '  .  Tfi  ^  +  iY2(u)^ 

(cos  cp) 


(83) 


.?{l  +  .09[to(2.5  C0h)  +  to  5T^pj-fi} 


(84) 


L.35^1  +  .07  Ui  ^  ^ ^ )  for  sunspot  minimum 


L.45^1  +  .06  in  2  5  %|  g  flj  ^or  BUn8Pot 
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Equations  (82)  and  (83)  have  the  same  form  as  Eq.  (80) .  Further¬ 
more,  since  2.5  CQ  a*»d  1*1  a*e  nearly  equal,  the  major  difference 

between  the  transmission  matrices  for  the  two  forms  of  v, /V^  lies  in 

i  a 

the  difference  in  the  numerical  value  of  V.  This  parameter  is  a  mea¬ 
sure  of  the  ratio  of  the  effective  size  of  the  transition  region  to 

that  of  the  Hall  region.  The  significance  of  its  magnitude  can  be  seen 

2 

.  To  compare  the  two  forms  of  v./V  ,  we  shall  con- 
x,y'  v  i  a’ 


by  examining  |T 


aider  cos  cp 

(80)  by  the  essentially  equivalent  parameter  2.5  uih£ 


1  (high  magnetic  latitude),  and  replace  1.1  ujhz^  of  Eq. 


The  magnitudes 


of  T 


x,v 


in  both  cases  can  then  be  written  as 


For  positive  u>  (transmitted  p-mode) ,  |t  j  has  a  maximum  whose 

y 

frequency  and  shape  depend  upon  Y.  For  small  Y>  the  peak  is  high  and 
narrow,  and  gives  rise  to  resonance  phenomena.  For  larger  values  of 
Y,  the  frequency  of  the  maximum  is  shifted  down,  and  the  peak  is  broad¬ 
ened  and  depressed.  We  can  see  from  Eqs.  (81)  and  (84)  that  for  the 
same  set  of  parameters  (C  ,  h,  U), ,  V  ) ,  the  long-range,  inverse- square 

2  O  13 

form  of  v,/V  glvos  values  of  Y  about  twice  aj  large  as  those  obtained 

1  9 

from  the  short-range  exponential  form.  Qualitatively,  of  course,  this 
could  be  predicted;  extending  collisions  to  higher  altitudes,  thereby 
increasing  the  size  of  the  transition  region  without  altering  the 
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effective  size  of  the  Hall  region,  will  increase  the  effect  of  collision 
broadening  of  the  transmitted  spectra. 

2 

In  Figs.  2  and  3  we  show  |t  (cp  -  G)|  for  both  forma  of  v. /V  . 

x,  y  i  a 

The  curves  for  the  exponential  form  differ  from  those  in  I  because  they 
are  based  on  a  different  set  of  parameters.  The  values  of  Y  used  In  I, 
about  .3  and  .45,  are  smaller  than  tho^e  derived  from  current  ionospheric 
profiles . 

The  period  of  the  transmission  maximum  is  given  by 


T 

m 


In 

<D 

m 


*v»v  9 

K0a  +  /> 


(87) 


8  seconds  at  sunspot  minimum  1 
22  seconds  at  sunspot  maximum ) 

15  seconds  at  sunspot  minimum  j 

I  ^2 

40  seconds  at  sunspot  maximum J  Vfl 


The  numerical  values  of  the  periods  given  in  Eq.  (87)  are  those  for 

high  latitudes.  The  period  would  be  somewhat  longer  at  lower  latitudes 

for  the  same  values  of  h  and  C  because  of  the  decrease  in  tv  . 

o  i 

The  maximum  value  of  |Tx  is  only  about  20  percent  larger  thau 
its  value  for  u)  -»  0  for  the  inverse-square  form  of  v./V^.  Qualitatively 
similar  results  apply  to  the  exponential  form.  Thus  the  transmitted 
p-mode  is  slightly  enhanced  over  a  wide  range  of  low  frequencies  from 
0  to  2/Tm.  At  middle-high  latitudes,  this  corresponds  to  an  enhance¬ 
ment  for  frequencies  up  to  about  .1  ops  at  sunspot  minimum  and  up  to 
about  .04  cps  at  sunspot  maximum;  thereafter,  |T  |  decreases  with 

y 

increasing  frequency. 

When  cos  cp  4  1,  |t  |  and  |t  |  are  equal  for  the  exponential  form 

2  x  y 

of  v,/V  ,  and  differ  only  slightly  from  each  other  for  the  inverse- 
square  form.  From  Eqs.  (82)  and  (83),  we  can  see  that  the  inclusion 
of  cos  cp  effectively  increases  C  (the  increase  of  T  in  the  inverse- 


Fig.  2  —  A  comparison  of  the  magnitudes  of  the  elemer'3  of  the  daytime 

magnetic  transmission  matrix  at  sunspot  minimum  for  two  forms  of 

ui /Va>  a  long-range  inverse  square  form  and  a  short-range 

exponential  form.  The  earth's  field  is  assumed  to  be 

nearly  vertical  (cos  cp  =“  1),  so  that  |l  |  |l  |. 

x  y 
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Fig.  3  —  A  comparison  of  the  magnitudes  of  the  elements  of  the  daytime 

magnetic  transmission  matrix  at  sunspot  maximum  for  two  forms  of 

Ui/vJ,  a  long-range  inverse  square  form  and  a  short-range 

exponential  form.  The  earth's  field  is  assumed  to  be 

nearly  vertical  (cos  9  “*  1),  so  that  | T  |  °*  | T  |  • 

x  y 
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square  case  being  larger  than  that  of  T^) .  Increasing  Q  decreases  the 
magnitude  of  |l  |  and  shifts  the  frequency  of  the  maximum  downwards. 

*>  y 

So  far,  we  have  discussed  only  the  transmitted  p-mode.  The  trans¬ 
mitted  e-mode  can  be  obtained  from  Eqs.  (80)  and  (81)  for  the  exponen- 
2 

tial  form  of  v. /V  ,  and  from  Eqs.  (80)  to  (84)  for  the  inverse  square 

1  O 

form,  by  replacing  iu  by  -iu  and  taking  the  complex  conjugate  of  T 

y 

The  magnitude  of  T  for  negative  u>  is  a  monotonically  decreasing 
x,  y 

function  of  |  w| . 

We  shall  postpone  a  discussion  of  the  transmitted  fields  for  an  in¬ 
cident  F-  or  S-wave  for  different  values  of  cp  until  we  have  found  solu¬ 
tions  for  the  higher  frequencies. 

The  approximation  for  small  <u  is  valid  for  frequencies  such  that 


a 


At  high  latitudes  (<u^  ~  200  rad/sec)  the  approximation  is  valid  for 

the  inverse  square  form  of  v.  for  frequencies  less  than  about  .23 

1  8 

cps  at  sunspot  minimum  and  less  than  about  .1  cps  at  sunspot  maximum, 
and  for  the  exponential  form  for  frequencies  less  than  about  1  cps  at 
sunspot  minimum  and  less  than  about  .25  cps  at  sunspot  maximum. 

SOLUTIONS  FOR  HIGHER  FREQUENCIES 

A  simple  means  of  estimating  the  transmitted  fields  at  higher 
frequencies  is  the  WKB  method.  This  approximation  is  applicable  to 
an  equation  such  as  Eq.  (69)  when  the  fractional  change  in  the  magni¬ 
tude  of  the  index  of  refraction  is  small  in  a  wavelength.  The  cri¬ 
terion  for  the  applicability  of  the  approximation  is  given  in  terms 
of  q^  by 
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1  iia  JLi  <  i 

2  'dZ  2'  1 

q 


(88) 


In  the  regions  in  which  this  condition  is  satisfied,  the  solution  of 
Eq.  (69)  can  be  approximated  by  a  linear  combination  of  the  functions 
t±(2),  where 

±i  J^qdz 

U*>  -  - -  (89) 

*  Vq 


2 

For  the  function  q  given  by  Eqs.  (64)  and  (65),  Eq.  (88)  is  satisfied 
if 


|4Vj£J  >  1 


1 8  u£o|  >  1  for 


a 


These  conditions  essentially  cover  the  frequency  range  in  which  the  low- 
frequency  approximation  breaks  down. 

The  application  of  the  WKB  method  -o  v  .^*ution  of  Eq.  (69)  in¬ 
volves  some  subtleties  (especially  in  the  case  of  negative  frequencies), 
which  are  discussed  in  the  Appendix.  In  this  section  we  present  the 
results. 

For  positive  frequencies  (denoted  by  the  subscript  p,  since  they 
correspond  to  the  propagating  mode),  the  approximate  solution  of  Eq. 

(69)  is  giver  by 


2 

a  -  cos  (p 

y 


(98) 


1 
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As  before,  the  choice  of  the  subscript  x  or  y  depends  upon  the  boundary 
conditions  on  the  incident  wave. 

At  the  bottom  of  the  ionosphere  (z  =  zp) ,  Eq.  (90)  becomes 


e 

P 


exp 


qp(u)du 


(99) 


From  Eq.  (99)  we  obtain  the  magnitude  of  the  positive- frequency  elements 
of  the  transmission  matrix: 


iTx(cp)lP  "  (cos  <p) 1/4  e*p[-  Ju  qpx(u)duJ  (100) 


|Ty<Cp)  1  p 


"jfa-  "  (cos  «p)3/4  exp|-  Im  Ju  qpy(u)duj  (101) 


where  Im  Jqdu  denotes  the  imaginary  part  of  the  integral. 

The  effects  of  collisional  absorption  are  contained  in  the  expo- 

1/4 

nential  terms  in  Eqs.  (100)  and  (101).  The  factor  U)  in  6(zp)  (which 

leads  to  the  factor  ®-3/4  In  Sground  an<*  *n  l^x  y^  1®  a  consequence 

of  the  increase  in  the  index  of  refraction  (for  uS  <  1)  by  a  factor  of 
•1/2  # 

u>  between  the  Alfven  and  Hall  regions.  The  amplitude  of  the  elec¬ 
tric  field  in  each  region  is  inversely  proportional  to  the  square  root 
of  the  region's  index  of  refraction. 

Evaluating  the  integrals  in  Eqs.  (100)  and  (101)  for  both  the  ex- 

2 

ponential  and  inverse  square  forms  of  v./V  ,  we  obtain 

X  A 
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IV^lp "  1374T  (cos  9)  6 


_pr  .  -  ■?  Q  \ (1  -  V a)  cos  9) 

v2  ,  .1/4  2  *0  » cos  cp  Y 

■■  ii  1  rr»o  rni  p  ” 


- <^/  T  « 

a)  n 


I  Ty(cP)  I  p  "  _'3/4  _  (cos  <P)3^4  e 


.  JL.  f  (1  .J-JL- ) 

2  "  cob  9  o\  "cos  9  / 


«>  h 


(102) 


^4  */C, 


(103) 


v'5’  -  VS 


iTx(<p)  I  W  --jy4-  -  (cos  cp)A/  e 


^  Lin  ~ 


at  h 


VT 


lT  («P>  I  »»  r374  -  Ceos  9)J/  e 


-  V2S  — - —  in  — ^ 
3/4  _  cos  9 


(104) 


(‘ '  t) 


(105) 


w'  h 


To  compare  the  results  for  the  two  forms  of  v^/Va>  we  let  cos  9~  1  and 
use  a  value  of  m  200  rad/sec  appropriate  to  high  latitudes. 

We  then  obtain 


sur  jpot 

|t  (o,r>  „-^e"7Vfd-  -»7V5 

*»y  p  f3/4 


sunspot 

!*,.,<•> r  -774  ''i-wr<i  •  •m/r) 


(106) 


V1  z/^o 
"2  *  e 


(107) 


(0)| 


sunspot 
1  min 


lTx,y(°)l 


sunspot 

max 

P 


-1.2  >/F(l  +  4 

e  5 


-2.4v^(l  +  \ 

e 


y± 

»2 


(108) 


(109) 


where  f  Is  the  frequency  In  cps . 
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For  both  forms  of  v, /V  ,  the  major  frequency-dependence  of  |lj 

1  a,  3/4  11 

is  provided  by  the  factor  1/f  .  This  is  due  to  the  fact  that  the 

exponential  term  in  |t|^  does  not  vary  much  in  this  frequency  range, 

indicating  that  the  collisional  attenuation  apparently  does  not  depend 

much  upon  frequency. 

Although  the  exponential  term  in  the  elements  of  the  transmission 

matrix  does  not  vary  much  with  frequency  for  either  the  exponential  or 

2 

inverse- square  form  of  v,/V  ,  its  value  is  different  in  the  two  cases, 

1.  fl 

the  attenuation  being  predictably  larger  for  the  long-range  than  for 
the  short-range  form.  The  fractional  difference  between  the  elements 
of  the  transmission  matrix  for  the  exponential  and  inverse-square  forms 
is  appreciable  (especially  at  sunspot  maximum)  at  frequencies  below  a 
few  cps,  but  becomes  smaller  as  the  frequency  approaches  the  ion  cyclo 
tr on- frequency.  This  is  a  consequence  of  the  "shrinking"  of  the  transi¬ 
tion  region  as  the  frequency  increases. 

When  cos  cp  4  1,  we  can  see  from  Eqs.  (102)  to  (105)  that 
| T  |  <  | T  |  The  ratio  of  |T  |  to  |T  |  is  smaller  for  the  inverse- 

yp  xp  yp  xp 

square  form  than  for  the  exponential  form,  and  for  the  inverse- square 
form  is 


l  •£  |  - 

w  (co8  <p) 1/2  |l  -  1.2  vT'l  +  |  l n  -^(sec  <p  -  l)j  (110) 


at  sunspot  minimum,  and 


If  1 

pyE  ca  (cos  cp)^|l  -  2.4>/f^l  +  fn  ~yj(sec  cp  -  l)j 


(111) 


at  sunspot  maximum. 

When  cp  ~  30  deg  and  f  ~  1  cps,  the  ratio  is  about  0.7  at  sunspo- 
minimum  and  about  0.5  at  sunspot  maximum,  and  decreases  as  cp  increases. 

For  negative  frequencies  (denoted  by  the  subscript  e,  since  they 
correspond  to  the  evanescent  mode),  the  approximate  solution  of  Eq. 

(69)  is  given  by 
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ee(v>  «  -  f  <BX  -  iBy  cos  9)lnc  exp 


KOw1  q*<v°,dv  ■ «)] 


X  cosjj 


*V  ■  qe(v)i  tt 

— — -  q  (v  )dv  -  7  I  for  v  »  v0 

v0  *e(vo>  c  °  41  2 


(112) 


£  (v)  w  (B  -  IS  cos  cp).  G)^4 

ev  '  c  x  y  T/lnc 


,2  +  £28jf| 


1/4 


a 

x.yl 


f  /  p  o  q  (v)  /2  q  (v)  \i 

“‘P[1(IV  “'<v°)dv '  ’e(v  %  '^V|dv)l  <U3) 


2.  q„(v), 


for  v  y  V2 

where 


Ill 

c 

* 

(114) 

ee(V)  = 

(115) 

qe(v)  =  q*(-«,u*) 

(116) 

and  is  defined  by 

qe(v2)  -  0  (117) 

From  ^e(v^)  we  cin  find  CB*(~«>)  ]Grou.nd  an<*  C^e  ne8a*ive- frequency  ele¬ 
ments  of  the  transmission  matrix.  The  results  9re  given  by 
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|ix«p)lp 


(U8) 


*i  z/Co 


|Ty<cp)  |  e  -  e 


_J_l£L  i 

”  cos  cp  1 


lTyWlp 


(119) 


T  (<p)| 

x't/  i  e 


-  fv^Visj  c. 


|t  | 

1  x'p 


(120) 
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a  1 1 


(‘-ir) 


|TyWle 


Sv^Viii  — — 

2  ^ 1  1  cos  cp  i_  | 

y  p 


(121) 


The  transmitted  amplitude  of  the  e-mode  falls  off  exponentially  with 

ij  | uS|  for  |(u|  larger  than  about  4  rad/sec  at  sunspot  minimum  and  for 

|(u|  larger  than  about  1  rad/sec  at  sunspot  maximum.  This  applies  to 
2 

both  forms  of  v. /V  ,  although  the  magnitudes  of  the  transmitted  ampll- 

X  fl 

tudes  differ  cot  the  two  cases  because  of  differing  amounts  of  colli- 
slonal  abc  -ption.  The  exponential  decrease  of  amplitude  with  y[\*\ 
is  a  consequence  of  the  evanescence  in  the  Hall  region,  and  is  not  the 
result  of  true  absorption.  In  Figs.  2  and  3  we  show  |t  |  when 

2  x»y  p>® 

cos  cp  ^  1  for  both  forms  of  v ,  /v. 

1  a 

Cur  analytical  approximations  (using  the  inverse- square  form  of 

2 

v^/Va)  are  compared  with  the  numerical  results  of  III  in  Figs.  4  and  5. 

They  agree  well  for  | T  |  at  frequencies  up  to  sbout  5  cps  and  for 

x  p 

jTyjp  at  frequencies  up  to  about  1  cps.  Our  approximations  give  larger 

values  for  |t  |  than  those  of  III,  the  error  becoming  fairly 

x»y  e 

large  et  the  higher  frequencies;  this  is  not  surprising,  in  view  of 
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g.  4  —  A  comparison  between  the  analytic  approximations  to  the 
elements  of  the  daytime  magnetic  transmission  matrix  obtained 
in  this  paper  (using  the  inverse  square  form  of  Vj/Vg)  and 
the  numerical  values  obtained  by  Karplus,  Francis ,  and 
Dragt  (1962).  The  magnetic  dip  angle  is  assumed  to 
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the  exponential  nature  of  the  frequency-dependence  of  | T | e .  However, 
except  under  certain  conditions,  which  we  shall  discuss  shortly,  the 
transmitted  field  consists  almost  entirely  of  the  p-mode  at  frequencies 
rbove  about  .5  cps,  and  we  need  not  be  concerned  with  the  small  amount 
of  e-mode  that  is  transmitted. 

Having  found  the  transmission  matrix  elements,  we  can  now  find 
|B|groun<j  for  an  incident  F-  or  S-wave  from  Eqs.  (48)  to  (57): 


Incident  F-Wave 


lBloro™<;-"/rlT<0>lp“<P-0  <124) 


'Ground 


-  VT  |T(0)  |  e  as  <p 


(127) 


In  Fig.  6,  we  show  l®l(;roun(j  aa  a  function  of  cp  for  incident  S- 
and  F-waves  of  frequencies  .5,  1,  and  2  cps  for  minimum  sunspot  con¬ 
ditions.  For  T  ,(cp)  in  Eqs.  (122)  to  (127)  we  have  used  the  values 

x>  y  , 

obtained  with  the  inverse  square  form  of  v^/V“. 
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VII.  DEPENDENCE  OF  DAYTIME  TRANSMISSION  CHARACTERISTICS 
ON  IONOSPHERIC  ELECTRON  DENSITY 


So  far,  we  have  expressed  the  elements  of  the  transmission  matrix 
In  terms  of  dimensionless  variables  whose  units  of  frequency  and  length 
are  given  by: 

u)io  (unit  of  frequency)  210  cos  <f  rad/sec  (at  nonequatorlal  latitudes) 

(128) 


V 

on 

-  (unit  of  length) 

10 


1.2 

vl/2 

^io5 


km 


(129) 


where 


Nio  -  N1q  in  (No./cc)  x  10" 6  (130) 

and  the  subscript  o  refers  toi>0  (the  altitude  at  which  -  Wj) . 

The  order  of  magnitude  of  Nlo  la  .1  under  normal  daytime  conditions. 

(In  our  calculations  ve  have  used  values  r*  N^q  of  .1  at  sunapot  mini¬ 
mum  and  .175  at  sunspot  maximum.) 

Let  us  examine  the  dependence  of  the  daytime  transmission  charac¬ 
teristics  on  N^q.  First,  we  express  the  quantity  *VJsic0  in  terms  of 
N.  ,  since  this  quantity  determines  which  of  the  approximate  solutions 

aO  mm 

cf  the  propagation  equations  is  applicable.  The  quantity  4^1  °o| C 
given  by 


(ui) 

where  is  in  km.  According  to  the  profiles  of  Sims  and  Bostick, 

appears  to  increase  with  N^q.  Using  the  numerical  values  for  minimum 

and  maximum  sunapot  conditions,  we  can  obtain  an  approximate  empirical 

relationship  between  C  and  N,  : 

o  io 
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fo 

10 


(1°  Hl0) 


3/4 


(132) 


which  gives 


Wiiic.  - 2  <10  »io)5/4  (U3) 

/_  _  5/4 

When  2vf(10  NiQ)  <  1,  we  use  the  results  of  the  low-frequency 
approximation.  We  have  seen  that  at  these  low  frequencies,  y  end 
h£Q/ui  are  the  quantities  determining  the  transmission  characteristics, 
particularly  the  period,  the  height,  and  the  width  of  the  p-mode  trans¬ 
mission  maximum.  Since  y  is  the  dimensionless  ratio  of  the  effective 
size  of  the  transition  region  to  the  size  of  the  Hall  region,  its  value 
is  independent  of  the  unit  of  length. 

For  T  ,  the  period  at  the  transmission  maximum,  we  have  obtained 
m 

[Eq.  (87)3  from  the  low-frequency  approximation 
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pa 


2nhCo  2.5(1  +  Y  ) 
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2*o 

COS  Cp 


(1  +  ■v2)»lo 


seconds 


(134) 


Since  the  transmitted  p-mode  is  enhanced  for  frequencies  up  to  about 

2/t  ,  larger  particle  densities  result  in  a  smaller  range  of  frequen- 
m 

cies  of  enhanced  transmission. 

r-  -  5/4 

When  2vf(10  N^q)  >  1,  the  positive  frequency  transmission 
matrix  elements  are  given  by 


774  * 


exp I-  1.2 


r~ 

f 

/  a.. 

cos  cp 

5l5<10Mio>1/2(1  +  f  toViiR]  <l35> 
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2 

The  exponent  In  Eq.  (135)  Is  based  on  the  inverse-square  form  of  v,/V  , 

1  fl 

which  seems  to  give  a  reasonably  good  fit  to  current  ionospheric  pro¬ 
files.  Since  the  exponent  in  Eq.  (135)  is  of  the  order  of  unity  for 
f  near  1,  small  variations  of  10  about  unity  will  not  have  much 

effect  on  I T I  . 

i  ip 

The  transmitted  e-mode  is  more  critically  affected  by  changes  in 

N.  .  For  It  |  ,  we  have 
io  1  x,y'e 


|T  |  M  exp  -1.1  J —  (10  N,  )1/2|t  | 

1  x,y'e  scos  cp  10  '  io'  1  x,y'p 


f3/4  .  « 

(io  N4  )a/  |t  for  f  near  1. 

.7  '  io'  1  r.y'p 


(136) 


Thus  the  attenuation  of  the  e-mode  increases  roughly  as  the  square  of 
the  attenuation  of  the  p-mode  as  increases.  When  cp  >  v/icu  (or  for 
dip  angles  less  than  about  75°  for  f  ~  1  cps) ,  the  e-mode  does  not 
contribute  much  to  the  transmitted  fields  and  we  need  not  be  concerned 
with  | T 1  e .  For  the  daylight  polar  regions,  appreciable  Increases  in 
Nio  would  result  in  a  sizeable  reduction  in  the  transmitted  amplitude 
of  an  incident  S-wave,  since  such  a  wave  gives  rise  mainly  to  a  trans¬ 
mitted  e-He. 

In  Fig.  7,  we  show  | B| Qroun<j  f°r  8  -spot  minimum  and  maximum  day¬ 
time  conditions  for  f  ■  .25,  .5,  1,  and  2  cps,  and  for  a  dip  angle  of 


74 


From  Fig.  7,  we  can  see  that  the  equation 


B 


Ground 


1.3 


e- 1.85 (10  Nlo) 


(137) 


gives  a  good  fit  to  the  theoretically  calculated  points.  We  can  use 
Eq.  (137)  to  estimate  | B | Qrountj  at  middle  latitudes  for  frequencies  in 
the  range  1/4  £  f  £  2  cps  and  for  values  cZ  not  much  less  than  0.1 
(sunspot  minimum)  or  much  greater  than  .175  (sunspot  maximum).  (Ex¬ 
trapolation  to  values  of  much  outside  the  range  shown  in  Fig.  7 


0.05  0.1  0.15  0.2 

Njo  (ion  dens!ly  at  z  =  0  in  (no/cc)  X10  6) 


Fig  7  --  The  amplitude  o£  the  daytime  ground-level  magnetic  signal  as  a  func¬ 
tion  of  the  ion  density  at  z  ■  0.  The  data  points  arc  obtained  from  the 
analytic  approximations  to  I B iGround.  derived  in  the  text;  the  nu¬ 
merical  values  of  the  parameters  are  from  Sims  and  Bostick 
(1963)  for  minimum  and  maximum  sunspot  conditions  and  a 
magnetic  dip  angle  of  74°  •  The  solid  curves  are 
empirical  fits  to  the  data  points  given  by 
Eq.  (137). 
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cannot  be  justified  without  more  information  on  the  variation  of 
with  N^q.) 

From  Eq.  (137)  we  see  that  for  a  given  frequency  in  the  frequency 
range  ci  .25  to  2  cps,  an  increase  in  electron  density  of  the  lower 
ionosphere  of  about  55  percent  results  in  a  reduction  of  the  transmitted 
amplitude  of  about  1/e.  The  shape  of  the  amplitude- frequency  curve  is 
Insensitive  to  because  of  the  insensitivity  of  the  colllsional 
attenuation  to  frequency.  This  implies  that  variations  in  daytire 
electron  density  (with  season,  time  of  day,  and  solar  activity)  should 
not  materially  affect  the  relative  transmitted  amplitudes  for  different 
frequencies  between  about  .25  and  2  cps,  unless  the  relative  amplitudes 
of  the  incident  fields  in  the  lower  Alfven  region  are  affected. 

Let  us  compare  Eq.  (137)  with  an  empirical  equation  obtained  by 
Wentworth  (1964)  from  the  numerical  results  of  III.  Wentworth  related 
the  attenuation  1/ 1 B | Qroun<j  to  t*ie  frequency  and  to  the  F2  ion-density 
maximum  by  the  equation 

1  £(foF2>4 

attenuation  2  rri -  -  1  +  — r^n -  (138) 

I "'Ground 


where  f  Fj  is  the  plasma  frequency  at  the  F2  ion-density  maximum  in 
megacycles/ sec,  and  is  given  by 


1.24  X  10 


4  <W‘ 


cm 


ion  density  at  F2  peak 


(139) 


When  f  >  .25  cps  and  F2  peak  ion  densities  are  greater  than  about 
.5  X  10^  per  cm^,  Eq.  (138)  becomes 
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Ground 


400 

*<w4 


(140) 


-52- 


The  frequency  dependence  of  Eq.  (140)  agrees  with  Eq.  (137).  A  further 
comparison  between  Eqs.  (137)  and  (140)  is  difficult,  since  Eq.  (  40) 
relates  l®iGroun(j  directly  to  the  F2  peak  ion  density,  whereas  Eq.  (1371 
relates  l®lgroun<j  to  t*ie  *-on  density  in  the  lower  ionosphere.  Since 
the  lower  ionosphere  is  the  region  in  which  most  of  the  collisional 
attenuation  occurs,  it  seems  preferable  to  us  to  relate  the  attenuation 
directly  to  ion  densities  in  this  region,  rather  than  to  those  at  the 
F2  maximum,  even  though  more  information  is  usually  available  concern¬ 
ing  the  latter. 

Our  theoretical  results  are  supported  by  recently  reported  exper¬ 
imental  observations.  Dawson  (1967)  reports  no  particular  correlation 
between  amplitudes  (of  disturbances  in  the  pc  1  frequency  range)  and 
ionospheric  parameters  at  College,  Alaska  (65°N) ,  in  agreement  with  our 
prediction  that  moderate  fluctuations  in  have  little  effect  on 

.5,  .75,  and  1  cps  made  at  Bale  St.  Paul,  Quebec  (47°N) ,  Campbell  (1967) 
reports  a  decrease  in  daytime  amplitudes  by  a  factor  between  1.5  and  3 
as  the  Xp  index  increases  from  0  to  The  relative  amplitudes  at  the 

different  frequencies  vary  little,  however,  with  the  Kp  index.  Further¬ 
more,  there  is  no  discernible  dependence  of  the  daytime  amplitude  ra¬ 
tios  on  the  solar  zenith  angle. 

So  far  we  have  considered  only  daytime  transmission.  In  the  next 
section  we  shall  consider  the  day--night  ratio  of  transmitted  amplitudes. 

i 


•Ground' 


From  observations  of  hydromagnetic  emissions  of  frequencies 
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VIII.  DAY-NIGHT  AMPLITUDE  RATIOS 

For  w^.ves  of  frequencies  below  c  few  cps,  the  "bottom?1  of  the 
nighttime  ionosphere  is  located  in  the  Alfven  region,  at  an  altitude 
of  about  300  km.  At  lower  altitudes,  the  Alfven  speed  is  sufficiently 
large  that  the  medium  behaves  essentially  like  a  vacuum- 

Within  the  nighttime  ionosphere,  the  fields  have  the  form  of  Eqs. 
(37,i  to  (AO).  If  the  origin  of  the  coordinate  system  is  taken  at  the 
base  of  the  ionosphere,  the  boundary  condition  on  E  at  z  ■  0  it  given 
by 


/dEx  y' 
\  dz  ) 


z-0 


"  ■  :  E,  Jz  ■  0) 


H 


x,y 


(141) 


where  H  is  the  height  of  the  nighttime  earth-- Ionosphere  cavity.  We 

can  then  obtain  E  (z  ■  0)  in  terms  of  A  and  A  : 

x,y  x  y 


Ex(z  -  0) 


cos 

2ik  H 
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Z 

ik  H  -  1. 

(142) 


Ey(2-0)  --~ 


2iC  H 
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Ilk  H  -  lJ 

x 


(143) 


Applying  the  induction  equation  at  the  surface  of  the  earth,  and  using 
the  boundary  conditions  [Eqs.  (48)  to  (52)]  to  give  A  and  A  ,  we  get 

x  y 


for  B 


Ground 


|  night 
'  Ground 


1  +  E2  H2 

*»y 


(144) 


If  we  use  the  parameters  appropriate  to  high  latitudes  obtained  from 
Sims  and  Bostick,  Eq.  (144)  becomes 
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B1  Ground  w  f  ^10  \o^ 


-1/2 


1  + 


4f2(10  Nio) 


1/2 


(145) 


The  approximate  day --night  amplitude  ratio  when  1/4  <  f  &  2  cps, 
0.5  g  ,«.*  .22,  and  tp  >  V5uj  is  given  by* 


I»IX 


-  1.3  exp  [-1.85(10  Hlo)][l0  Hl0]1/2 


1  * 


4? (10  Nlo) 


li 


(146) 


(We  have  let  coa  cp  -  1,  since  cos  cp  varies  only  slightly  for  values  of 
cp  appropriate  to  middle  and  high  latitudes,  and  the  dependence  on  cp  is 
furthermore  obscured  by  the  uncertainty  in  the  value  of  H.)  Let  us 
note  that  daytime  collisiunal  absorption  is  not  the  only  factor  influ¬ 
encing  this  ratio;  the  ratio  of  the  heights  of  the  effective  nighttime 
and  daytime  earth- -ionosphere  cavities  also  plays  an  important  role. 

We  can  see  from  Eq.  (146)  that  for  frequencies  larger  than  about 

—  1/2 

1/2(10  N^q)  ,  the  day — night  amplitude  ratio  becomes  almost  indepen¬ 
dent  of  frequency.  Figure  S  shows  the  day--night  amplitude  ratio  for 
frequencies  in  the  range  .25  £  f  £  2  cps  for  values  of  of  .05,  .1, 
.15,  and  .2. 

These  results  can  be  compared  with  the  experimental  observations 
reported  by  Ccmpbell  (1967)  on  the  relative  signal  amplitudes  of  hydro- 
magnetic  emissions  received  at  magnetically  conjugate  stations  during 
periods  of  changing  ionospheric  illumination.  The  stations  studied 
ere  at  Eights,  Antarctica  (75°S,  77°W)  and  Bale  St.  Paul,  Quebec  (47°N, 
70°W).  Events  occurring  simultaneously  at  the  two  sites  were  recorded, 
with  a  preferred  selection  of  events  which  occurred  vhen  the  Antarctic 
station  refined  in  darkness  as  its  Canadian  conjugate  progressed 
through  day— night  cycles.  The  day — night  amplitude  ratios  obtained 
ranged  from  .2  to  .7,  depending  upon  the  Kp  index  for  low  Kp  indices 
between  0  and  and  for  frequencies  between  .5  and  1  cps.  These 

In  Eqs.  (145)  and  (146),  refers  to  its  daytime  value. 


Fig.  8  --  The  ratio  of  the  amplitude  of  the  daytime 
to  the  amplitude  of  the  nighttime  ground-level 
magnetic  signal  for  a  given  incident 
field  amplitude  as  &  function  of 
frequency  and  daytime  ion 
density  at  z  ■  0. 


-56- 


results  are  entirely  consistent  with  our  theoretical  results  for 
.5  <;  1C  N^o  <1.5,  especially  considering  the  uncertainty  of  the  parame¬ 
ters  used  In  the  theoretical  calculations.  Furthermore,  the  pulsations 
measured  had  passed  through  the  lower  exosphere  and  the  F^  ion-density 
peak,  and  the  filtering  action  of  this  region  differs  somewhat  between 
night  and  day.  Finally,  our  calculations  are  based  on  the  assumption 
of  vertical  incidence,  whereas  the  observed  waves  apparently  traveled 
along  field  lines  between  conjugate  sites  and  were,  therefore,  obliquely 
incident  on  the  ionosphere.  One  can  expect  that  a  circularly  polarized 
Alfven  wave  propagating  along  the  field  lines  toward  a  middle- latitude 
ionosphere  will  break  up  into  an  S-wave  (E^  «  0)  and  an  P-wave  (B^  ««  0) 
iu  the  upper  ionosphere,  producing  a  mixture  of  the  two  types  incident 


should  have  a  value  somewhere  be- 


on  the  ionosphere.  If  so,  | B | Qroun<i 
tween  its  values  for  the  incident  S-  and  F-waves. 

In  any  case,  we  can  conclude  that  daytime  colllslonal  absorption 
does  not  cause  much  amplitude  attenuation  under  normal  conditions,  ac¬ 
counting  for  a  reduction  in  the  transmitted  amplitudes  by  a  factor  of 
about  3  at  f  ~  1  cps.  Daytime  *mplitudes  are  further  reduced  by  the 
transformation  of  part  of  the  incident  wave  into  the  highly  attenuated 
u-mode  in  the  lower  ionosphere.  However,  these  effects  are  partly  off¬ 
set  by  the  smaller  daytime  thickness  of  the  earth — ionosphere  cavity, 
which  increases  the  daytime  transmitted  amplitudes  relative  to  their 
nighttime  values.  The  outcome  of  these  competing  effects  is  a  day-- 
night  amplitude  ratio  of  the  order  of  .5  under  near-minimum  sunspot  con¬ 
ditions.  This  permits  a  projection  such  as  that  made  by  Campbell  (1967) 
of  ground  measurements  of  the  occurrence  of  pc  1  oscillations  and  of 
their  frequency  variations  into  the  upper  ionosphere  just  below  the  F£ 
peak  (at  nonpolar  latitudes).  It  is  not  possible  to  make  a  similar  pro¬ 
jection  of  ground-polarization  measurements  at  nonpolar  latitudes.  The 
fields  transmitted  through  the  daytime  ionosphere  should  be  predominantly 
p-mode  (right  circularly  polarized),  regardless  of  the  polarization  of 
the  incident  wave. 

Finally,  It  should  be  pointed  out  that  although  the  transmitted 
amplitude  should  not  undergo  large  diurnal  variations  (for  a  fixed  in¬ 
cident  amplitude),  the  fraction  of  the  incident  power  that  is  absorbed 
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may  differ  significantly  between  night  and  day.  The  fraction  of  the 
incident  power  ahsorbed  by  the  ionosphere  can  be  determined  from  the 
elements  of  the  reflection  matrix  (with  l  perfectly  reflecting  earth, 
the  energy  that  is  not  reflected  is  absorbed  by  the  ionosphere) . 

While  the  approximate  solutions  of  the  propagation  equations 
derived  in  this  paper  are  adequate  for  the  purpose  of  estimating  the 
transmission  matrix  elements,  more  accurate  solutions  are  needed  for 
the  determination  of  the  reflection  matrix  elements.  These  will  be 
derived  in  a  later  paper. 


BLANK  PAGE 


-59- 


APPEHDIX 


In  order  to  apply  the  WKB  method  to  the  solution  of  Eq.  (64),  it 

is  convenient  first  to  introduce  a  new  variable,  u,  in  terras  of  which 
2 

the  phase  of  q*"  (u)  is  constant,  except  at  possible  turning  points  where 

2 

its  s  gn  changes.  Tne  variable  u  and  the  function  q^Cu)  for  the  two 
forms  of  . /\)2  are  given  by 
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where 


a  =  1 
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a  =  cos  cp  (A-5) 
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We  consider  the  p  and  e  inodes  separately  for  both  forms  of  -v  /V^,  re¬ 
placing  Eq.  (69),  in  which  r  can  have  either  sign,  by  the  following  two 
equations  in  which  the  sign  of  jl  is  taken  to  be  positive; 
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(A- 10) 

The  function  q  (u)  is  given  by  (A-2)  or  (A -4),  depending  upon  the  func- 
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tional  form  of  v. /V  . 
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For  the  p-mode,  q^(u)  has  no  zeroes,  and  the  application  of  the  WKB 
method  to  the  solution  of  (A-6)  is  straightforward.  In  applying  the 
boundary  condition  Eq.  (34)  at  s  ,  we  neglect  1/Vouh  as  compared  to 
unity.  The  WKB  approximation  neglects  terms  of  order  1/VuiC  as  com¬ 
pared  to  unity,  and  since  h  ~  10  C  ,  it  would  be  inconsistent  to  re- 

o 

tain  terms  of  the  order  l/v/Sf!.  Thus,  Eq.  (34)  is  replaced  by  the 
condition  that  d£/dz  vanishes  at  z^.  This  gives  for  £  (u) 
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u0  “  u<*0)  (A- 14) 

U1  ’  u(*i> 

2 

For  the  e-mode,  qg(v)  has  a  zero  at  v^,  where 

1  z/CD 

v2  -  k i  ^  for  -  e  (A- 15) 


The  WKB  method  cannot  be  used  near  this  turning  point,  but  can 

be  applied  for  values  of  v  not  too  near  v„.  The  usual  technique  for 

L  2 

solving  this  t>pe  or  problem  is  to  replace  qg(v)  by  a  linear  function 
of  v  near  the  turning  point.  The  solution  of  (A- 7)  when  v  is  near  v2 
then  becomes  a  linear  combination  of  Bessel  functions  of  order  1/3: 

v0  ■  lciJi/3(u)  +  v-i/sH  (a-i7> 


where 
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When  |w|  is  large,  the  asymptotic  form  of  (A-17)  is  a  linear  combina¬ 
tion  of  WKB  solutions,  so  that  (A-17)  can  be  used  for  the  whole  range 


of  values  z  £  z  £  z, . 

o  1 


We  can  find  C ^  in  terms  of  by  imposing  the 


condition  that  (A-17)  should  give  a  decreasing  exponential  in  the  di¬ 
rection  away  from  the  turning  point  towards  the  region  of  evanescence. 
We  then  use  the  boundary  condition  at  vq  to  determine  C^.  The  applies* 
tion  of  this  procedure  is  straightforward,  and  gives 


ST 
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